Inelastic electron tunneling spectroscopy of a single nuclear spin 



O 
(N 



F. Delgado and J. Fernandez-Rossier 
Departamento de Fisica Aplicada, Universidad de Alicante, San Vicente del Raspeig, 03690 Spain 

(Dated: January 26, 2011) 

Detection of a single nuclear spin constitutes an outstanding problem in different fields of physics 
such as quantum computing or magnetic imaging. Here we show that the energy levels of a single 
nuclear spin can be measured by means of inelastic electron tunneling spectroscopy (lETS). We 
consider two different systems, a magnetic adatom probed with STM and a single Bi dopant in a 
Silicon nanotransistor. We find that the hyperfine coupling opens new transport channels which 
can be resolved at experimentally accessible temperatures. Our simulations evince that lETS yield 
information about the occupation of the nuclear spin states, paving the way towards transport- 
detected single nuclear spin resonance. 
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Probing a single nuclear spin represents both the ulti- 
mate resolution limit of magnetic resonance imaging and 
a requirement in quantum computing proposals where 
the nuclear spin is used as a qubit. The idea of stor- 
ing and manipulating information in nuclear spins goes 
back to quantum computing proposals based on P donors 
in SiQ and NMR quantum computing [2|. Because of 
their very small coupling to their environment, the nu- 
clear spin coherence time is expected to be very long but, 
for the same reason, quantum measurement of a single 
nuclear spin remains a formidable task. Recent experi- 
mental breakthroughs have made it possible to perform 
single shot non-destructive measurement of a single nu- 
clear spin by means of optically detected single spin mag- 
netic resonance in NV centers in diamond [3]. 

Here we propose a setup based on inelastic electron 
tunneling spectroscopy (lETS) that would permit to 
probe the spin transitions of a single nuclear spin with a 
space resolution down to 1 A outperforming in this par- 
ticular regard the optical detection. Our proposal builds 
on recent progress to probe the spin of a single atom using 
two different strategies. On one side. Scanning Tunneling 
Miscrocope (STM) inelastic electron tunnel spectroscopy 
[3-[loj which allow to measure the electron spin spectral 
function of a single magnetic atom [11] weakly coupled 
to a conducting substrate. On the other side, the fabri- 
cation of a Silicon nanotransitor where transport occurs 
through the electronic states of a single dopant (12l Il3|. 

In STM-IETS experiment electrons tunnel between the 
tip and the conducting substrate going through the mag- 
netic atom. As the bias voltage V is increased, a new 
conduction channel opens whenever eV is larger than the 
energy of some internal excitation of the atom. In the 
case of isolated transition metal atoms with partially full 
d shell, like Mn, Fe or Co, the only internal excitations 
available in the range of a few meV are spin excitations 
associated to the magnetocrystalline anisotropy [1, . 
In the case of the single dopant nanotransistor, the lETS 
of the electron spin in the donor level could be performed 
in the cotunneling regime [l3|- 



Exchange coupling to nearby magnetic atoms affects 
significantly the spin excitation spectrum of the atom 
under the tip 

S i, SEI- Hyperfine coupling to the 
nuclear spin should also result in a modification of the 
electronic spin spectral function which, in turn, could be 
probed in lETS provided that the spectral resolution is 
high enough. Our mechanism differs from earlier theory 
work 

proposing to detect the nuclear spin look- 
ing at its influence on the STM current noise spectrum. 
Their approach is based on previous experiments where 
electronic spin fluctuations are detected in the current 
noise spectra [l7|. not in the conductivity spectra. Prob- 
ing the spin transitions of a single nuclear spin would 
yield a completely unambiguous chemical identification 
of the atom and would be a first step towards transport- 
based quantum measurement of a single nuclear spin. 

The rest of this letter is organized as follows. We 
first describe the general theory that relates nuclear spin 
fiips to lETS transport features. Then we consider the 
archetypical case of a single Mn atom in a CU2N sur- 
face Q ■ We find that the detection of the nuclear 
spin excitations could be done at 4mK, below the recently 
demonstrated lOmK experimental limit [isj . and how the 
visibility can be enhanced driving the nuclear spin out of 
equilibrium. We then analyze the case of ^°^Bi in Sili- 
con and we find it is an optimal system to observe single 
nuclear spin fiips at temperatures up to 60mK. 

The electronic spin S and nuclear total angular mo- 
mentum / are described by a Hamiltonian 'Hq{S^ I) whose 
eigenvalues and eigenvectors are denoted by cm and 
1^) = E/,,s,*Af(/.,5,)|/,)|5,), where \h)\S,) is the 
basis in which both electronic and nuclear spin have well 
defined projection along the z axis. The spin mixing coef- 
ficients m{Iz, Sz) depend on the specifics of the Hamil- 
tonian, described below, which includes hyperfine cou- 
pling, Zeeman coupling and magnetic anisotropy terms. 
The inelastic transport spectroscopy is sensitive to tran- 
sitions between states M and M' with excitation energy 
^M',M = fA/' — Em- The electronic spin is coupled to 
electrons in the electrodes, denoted as tip (source) and 
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the surface (drain) in the STM (nanotransistor) geome- 
try, through a Kondo-hke Hamiltonianflll. [l9l-[2l| : 
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where both electrode conserving and electrode non- 
conserving exchange couplings are included. The op- 
erator c\ ^ creates an electron with spin a and orbital 

quantum number A — ri,k, where t] — T, S labels the 
electrode and k the wave vector. The index a runs over 
a = x,y, z, with t^"^ and Sa the Pauli matrices and the 
atom electronic spin operators respectively. The a = 
term (with r^''^ and 5*0 the identity matrix) corresponds 
to potential scattering. We assume that exchange is mo- 
mentum independent, spin isotropic and electrode de- 
pendent: Tx^y^a = VrjVjjiT- Here v^i are dimensionless 
parameters that account for the different coupling be- 
tween the magnetic adatom and either the tip or the 



surface 21, 22 



The inelastic current due to electronic-spin assisted 
tunneling /jn can be expressed [ll| as a convolution of 
the electronic spin spectral function: 
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where Pm denotes the average occupation of the state M 
and S^'^'^' = {M\Sa\M'). The inelastic current can be 
written as 
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Here go = ^GqPtPs \Tvtvs\ is the tunneling conduc- 
tance. Go = 2e^//i is the quantum of conductance and 
Pti the density of states of the rj electrode at the Fermi 
level. The inelastic current associated to a single channel 
is given by (e/Go)«(A + eV) = ^(A + eV) - ^(A - eV), 
with g{uj) = cj (1 - e-^'^y^ and 13 = l/keT. 

Importantly, '^^''^y'^^'' has a step at eV = A that 
accounts for the characteristic dl/dV lineshape. Thus, 
whenever the bias energy eV exceeds a transition energy 
between spin states with Am',m — eV such that the ini- 
tial state is occupied, Pm > 0, and the electronic spin flip 
transition is permitted, S^'^ ^ 0, the differential con- 
ductance dl/dV has a step and d^I/dV"^ has a peak (or 
a valley at negative bias) with a thermal broadening [23j 
of TksT = 5AkBT. Thus, the d^I/V^ line-shape should 
be quite similar to the electronic spin spectral function 
S{eV). In the following we assume that the intrinsic 
broadening of the spin excitations, due for instance to 
their coupling to the conducting substrate 21, 22, i^l, is 
negligible compared to Tk„T, as it happens in the case of 
Mn adatoms in CU2N [i,[2i,[2i|. 
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FIG. 1: (a) Energy level scheme of the ^^Mn^^ adatom on the 
CU2N surface, (b) Detail of the Sz = 5/2 — >■ 3/2 transition 
with resolved hyperfine structure, (c) dl* /dV and d^I* /dV'^ 
in the high current regime {vt = us = 1) at T = 4mK 
(J* — I/go). Inset: magnified low bias d? I* / dV^ . (d) Cor- 
responding occupations of each of the 12 lowest energy elec- 
tronuclear eigenstates of "^0(8,1) for the ^^Mn^^ versus ap- 
plied bias. 



The occupation functions Pm can differ substantially 
from those of equilibrium when the typical time elapsed 
between inelastic current events is shorter than the spin 
relaxation time [21], We determine the occupa- 



tion functions Pm by solving a master equation [2l|, [22 1 



that accounts for the dissipative dynamics of the current 
driven electronic spin interacting with the nuclear spin. 
Both energy and spin are exchanged between the elec- 
trons in the electrodes on one side, and the nuclear and 
electronic spin of the single atom on the other. The scat- 
tering events involve the creation or annihilation of an 
electron-hole either in the same electrode, in which case 
no current is involved, or in different electrodes. Thus, 
Pm depends in general, on the voltage, the conductance 
and the temperature 22 1. 

We now address the main question of this letter: un- 
der which conditions would lETS reveal transitions that 
provide information about the nuclear spin state?. We 
first consider the case of a single Mn adatom in a CU2N 
surface. This system has been widely studied experimen- 
tally and theoreticallyd i, i, HI Hfl-il [H [if but the 
role of the nuclear spin has been overlooked so far. The 
single atom electronic spin can be described by means of 
a spin S — 5/2 Hamiltonian, 



Hs = DSl + E{Sl - Si) + g.pBB ■ S 



(4) 



where D = — 39^eV and E — 7/ieV account for the uni- 
axial and in-plane anisotropy respectively, whereas the 
third term describes the electronic Zeeman coupling. At 
zero magnetic field the experimental 6] dl/dV features a 
step at an energy of A ~ 4|Z3| associated to electronic 
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spin flip between the two ground states which neglect- 
ing E, have Sz = ±5/2 and the first excited states, with 
^2 = ±3/2, see Fig.[T][a). The only stable nuclear isotope 
of Mn is ^^Mn and has a nuclear spin of / = 5/2, so that 
the electronic-nuclear system has 36 states in total. The 
hyperfine structure associated to the coupling between 
/ and S has been resolved with EPR and NMR experi- 
ments that address ensembles of more than 10^^ atoms, 
and can be described by a Heisenberg-type coupling. So, 
the spin Hamiltonian reads as 



■HoiS,l) = ns + AS.I 



(5) 



The strength of the hyperfine coupling A depends both 
on the nuclear magnetic moment and on the shape of the 
electronic cloud, which is environment dependent. In the 
case of Mn, A varies between 0.3 and 1/ieV [1^. Here 
we take A — IjieV . The effect of the hyperfine coupling 
is to split each of the 6 electronic levels into 6 nuclear 
branches, as observed in Fig.[T^,b). The lowest electronic 
multiplet, corresponds to the Sz = ±5/2, while the Sz — 
±3/2 and Sz = ±1/2 branches are found approximately 
at A\D\ and 6|D| above. 

The discussion can be simplified if we take advantage 
of the fact that \D\ ^ E^A so that the eigenstates of 
^0(8,1) have, to zeroth order in E and A, well defined 
projection of the electronic {Sz) and nuclear (Jz) spin, 
although the numerical calculations are done with the 
exact states. The energies are approximately given by 
Eo{Sz,Iz) = DS^^ + ASzIz + gt^sSz- We find two kind 
of transitions: low energy excitations with ASz — and 
higher energy excitations with AIz — 0,ASz — ±l.The 
former are found at small bias eV ~ SzAAIz, where AIz 
is the change in nuclear spin. Whereas these transitions 
yield very weak peaks in (Pl/dV^ , shown in Fig. Id), 
they contribute to drive the nuclear spin out of equilib- 
rium, as seen in Fig. 3a). 

In contrast, the electronic spin transitions between the 
electronic ground states ±5/2 and the first excited states 
±3/2 that, to zeroth order, conserve the nuclear spin, oc- 
cur at higher eV ~ 4|_D| and have a much stronger signal 
(Fig. lc,d) and 2). At i? = their energies are given by 
A±{Iz) = A\D\zLAIz so that, the hyperfine coupling splits 
the lowest energy line of the electronic spin spectral func- 
tion into 6 lines separated by ^ = A±{Iz) — A±{Iz ± 1). 
In order to have a thermal broadening smaller than A 
, fcsT must be reduced down to ~ 0.2^eV (T ~ 2mK), 
which would deplete the thermal occupation of the higher 
energy nuclear states within the electronic ground state, 
and the visibility of the corresponding spin excitations. 
In contrast, by decreasing the tip-atom distance [9] it is 
possible to drive the system out of equilibrium and pop- 
ulate also the higher energy nuclear spin states at low 
temperatures, as shown in Fig. Id), which makes it pos- 
sible to observe all the transitions, as seen in Fig. [Eleven 
at r = 4mK. 
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FIG. 2: (Color online) Conduction spectrum d?I* /dV'^ as 
a function of applied bias (/* — I/go), (a) Spectra for T — 
4mK (thick solid line) and T — 2mK (thin solid line) with 
vt = vs ~ 1 (far from thermal equilibrium) and T — 2mK 
(thick-dashed line) with IOut — vs = 1 (close to thermal 
equilibrium), (b) d^I/dV^ for T = 4mK and 
(solid line) and spin spectral function S{eV) as a function of 
applied bias (dashed line). 



In figure [^h) we plot d'^I/dV^ together with the the 
electronic spin spectral function S{eV), where the oc- 
cupations Pm evaluated at eV and the delta function 
replaced by ^ . It is apparent that d'^I/dV^ and S{eV) 
are related. Importantly, the spin spectral function con- 
tains information not only about the energy levels of the 
joint nuclear-electronic spin, but also about the occupa- 
tions of the states. Thus, the height of the d^I/dV^ can 
be correlated with the occupations of the nuclear spin 
states, outperforming current- noise spectroscopv fisi [l6| . 
This indicates that it would be possible to use lETS-STM 
as a detector in a magnetic resonance experiment, in 
analogy with the optically detected single spin magnetic 
resonance. This would permit to probe the magnetic field 
with subatomic spatial resolution and accuracy afforded 
by the linewidth of nuclear resonance experiment. 

Nuclear spins with larger hyperfine coupling, like 
Yb, Er or Pr, [13, ^ could be probed with lETS at 
higher t emp eratures than Mn. For instance i^'^Er'^"'" in 
¥28105 [23 can have A up to 6/ueV. The case of ^^^Bi in 



Silicon is particul arly interesting. This system attracts a 
lot of interest 2^ 301 in the context of quantum comput- 
ing based on the nuclear spin of donors [jy] where address- 
ing the nuclear spin of a single dopant is required. The 
hyperfine coupling between the 7 = 9/2 nuclear spin and 
the electronic spin of the donor state is quite large, A = 
G.l^ieV (1.48GHz)[2i Is^. The zero field Hamiltonian 
AI-S can be diagonalized in terms of the total angular op- 
erator F, resulting in two multiplets (F=4,F=5) with en- 
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FIG. 3: (Color online) Energy spectrum of ^"^^Bi in Silicon as 
a function of applied field and the corresponding spectra 
at ksT = WmK. 



ergy E{F) = 4 {F{F + 1) ~ S{S + 1) - 1(1 + 1)). The 
zero field splitting is as large as 5A ~ 30^eV, which 
could be resolved in lETS at fc^T ~ 60mK. The evo- 
lution of the spectrum as a function of B is shown in 
Fig. 3. A single Bi dopant in Si could be probed by single 
spin lETS-STM , recently demonstrated in semiconduc- 
tor substrates [lojl or in a single dopant Silicon nanotran- 
sistor both in the sequential regime, for which the 

dl/dV curve yield the single electron spectral function, 
with peaks at every one of the 20 energy levels [3l| , or in 
the cotunneling regime, described with the effective cou- 
pling (HI) which yields information of the electronic spin 
spectral function. In figure 3 we show the evolution of 
the (Pl/dV^ curves, as measured either with STM or in 
a nanotransistor in the cotunneling regime, for different 
values of B at fc^T" = lOmK. It is apparent that trans- 
port permits to probe the spin-flip transitions between 
the Zeeman-split states of the joint electron-nuclear spin 
system. 

In summary, we propose an experimental approach to 
probe a single nuclear spin using lETS of the hyper- 
fine structure of the electronic spin excitations, extending 
thereby the range of applicability of lETS. Our simula- 
tions show that this technique yields information both 
about energy levels and the occupation of the nuclear 
spin states of a ^^Mn adatom probed with a STM at 
4mK. In the case of a single ^^^Bi dopant in a Si nan- 
otransistor, the hyperfine structure could be detected at 
60mK, well within range of current state of the art. 
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